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Abstract. We present measurements of galaxy clustering at redshift (z) = 0.6 using 10 360 galaxies with photometric red- 
shifts over an area of 0.78 deg 2 from the COMBO-17 survey. To obtain a result that is unaffected by redshift uncertainties, we 
calculate the projected correlation function w(r p ), giving results for red sequence and blue cloud galaxies separately. The corre- 
lation function of the red galaxies displays clear deviations from a power law at comoving separations around 1 to 3 h~ x Mpc, 
and similar but weaker trends are suggested by the data for the blue galaxies. To interpret these results, we fit the correlation 
functions with analytical predictions derived from a simple halo occupation model. This combines linear clustering of the un- 
derlying mass with a description of the number of galaxies occupying each dark-matter halo (the halo occupation distribution). 
If the occupation numbers are taken to be a simple power law N oc M a , then a ~ 0.5 and a ~ 0.2 for red and blue galaxies 
respectively. These figures are little different from the values required to fit present-day clustering data. The power-spectrum 
shape is assumed to be known in this exercise, but we allow the data to determine the preferred value of as, the linear power- 
spectrum normalization. The average normalization inferred from red and blue galaxies at (z) = 0.6 is as = 1.02 ± 0.17 at 
zero redshift, consistent with independent estimates of this local value. This agreement can be regarded as a verification of the 
hierarchical growth of the halo mass function. 
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ISheth & Tormenl 19991k 

However, a long-standing challenge has been to under- 
stand how these ideas could be implemented in the context 



In current models of galaxy formation, structure grows hierar- 
chically from small Gaussian density fluctuations. Galaxies are 

presumed to form within virialized dark matter haloe s when ,. „ , , „ . . . , . , „. , , . 

, , , , t-— 71 of Cold Dark Matter (CDM) models. The galaxy correlation 

the baryomc gas cools and condenses into stars (e.g. Cole et al. . . , , , , , , , 

trr—t i_ ° , , . „ , . . TT~; , function has long been known t o be e xtremely close to a sin- 

2000). The formation and evolution of galaxies should thus be . . — .. . — || lr ,^nlln ill . 

, — , . , ... r, , , , ™ • gle power law (Totsuu & Kihara 1 969, Peebles 1 974), and yet 

closely tied to the merging history of dark matter haloes. This .. . , ', ,. ... 

_ , , . ,. this is not the case tor the nonlinear mass correlations in a 

paper uses measurements of galaxy clustering at intermediate , . . , , . . , 

, , ..... CDM model. Here, the matter correlation function rises above 

redshift to test this basic picture. . „ . . ^ . ,^„,, 

_ , . , , . , , , a best-fit power law on scales r & 1 h Mpc and falls below 

The complex relation between galaxies and dark matter . ^ . _ , _■, tz — r~r^ — , . II, ,,„ n | , 

,,,"..„ , • ,. it again on scales r 0.2 h Mpc (Jenkins et al. 1998 and 

has become clearer only slowly. Empirically, galaxies dis- ° . . . „. . , , , , . . , 

, , . . , . , , references therein). This puzzle was only resolved when it be- 

play biased clustering in which the amplitude of their cor- , , , • r , , , , 

, . . .-,1 iii. came clear that the correlation function of dark matter haloes 

relations varies with galaxy type: older galaxies are gen- .. , ,. ... .... . . , . . 

, , , , . (including subhaloes inside large host haloes) differs sigmfi- 

erally much more strongly clustered than young, starform- , „ , , . . . ^ , T . . 

, . 11 , cantly from the correlation function of the mass. In practice, the 

ing galaxies, and bright galaxies are more strongly clustered ,. . . . . . . . , , , „ „ 

. „. , . , r- — : — " — 1| n _ i i " ' -| predicted correlation function of galaxy-scale haloes follows a 

than faint galaxies (e.g. Davis & Geller 1976; Norberg et al. , . ,„„,_i, ., ... , . 



— — — , - ' IT ^^^L T ' — T, — 7r 7 power law down to 100 h 1 kpc, with an amplitude and slope 

2002; Phleps & Meisenheimer 2003). It is generally believed .. .... , , ■ — . r 7 v . — nr lnnn 

- — , — " — ■ ■ ■ — ~, , , , „ similar to the data on real galaxies (Kravtsov & Klvpin 1999; 

that such trends can be understood through the tendency for r-; : — ■ — — Tl U„„ .[ t^~~~ , , I L„„ .I U . . , Si 

,.ii i . i ^ J Nevnnck et al. 2004; Kravtsov et al. 2004; Tasitsiomi et al. 

dark-matter haloes to display clustering that is larger for rare L^J .1 ™ . . ... „ ., , . . 

2004). This phenomenon underlies the considerable scale- 
dependent bias predicted by semianalytic and hydrodynamic 
simulation models of galaxy formation rtColfn et alJ fl999: 
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These developments in turn stimulated a simpler and 
more direct insight into bias and its depende nce on scale, 
through the so-called halo model (e.g.LTingll998t ISeliakl2 000 ; 
IPeacock & Smithll2000t ICoorav & Shethl l2002 and references 
therein). Here, the shape of the correlation function is deter- 
mined by the linear clustering of the dark matter, and the rela- 
tion of the galaxies to the dark matter halos in which they reside 
(the Halo Occupation Distribution; HOD). In particular, a break 
in slope is expected when the correlation function changes from 
being dominated by pair counts of galaxies in separate dark 
matter halos to the small-scale regime, where pairs come from 
two galaxies that reside in the same halo. Any pure power law 
correlation function would require coincidental alignment of 
these two terms, and indeed analyses of the two-point correla- 
tion function of galaxies in the local universe have detec ted 
small deviations from the power-law form ( Hawkins et alJ 
l2003aUZehavi et alJ2004Eo05tlAbazaiian et all2005l) 

In this paper, we use the COMBO-17 survey dWolf etalJ 

2004) to carry out a similar investigation of the exact shape 
of the correlation function at higher redshifts. We calculate 
the projected correlation function w(r p ) for red sequence and 
blue cloud CO MBO-17 galaxies (following the definition of 
iBell et aljEooH) . in the redshift bin 0.4 < z < 0.8. By com- 
paring these results to the predictions of the halo model, we are 
able to infer the mean number of galaxies per halo of a given 
mass (the halo occupation number) and also the z — power- 
spectrum normalization as (the rms density variation averaged 
over 8ft." 1 Mpc spheres). 

This paper is structured as follows: The COMBO-17 survey 
and the data used in this analysis are briefly described in Sect. 
13 The halo model is introduced in Sect. [3] The method used to 
estimate the projected correlation function is explained in Sect. 
|4] In Sect.[5]we investigate the shape of the correlation function 
for red sequence and blue cloud galaxies, and in Sect.[6]the re- 
sults are discussed. We a ssume a cosmological geom etry taken 
from the WMAP results JSpergel et alJl2003ll20Q6h and the fi- 
nal 2dFGRS power spectrum results (Col e et alJl2005h : a flat 
model with f2 m = 0.25. All lengths quoted are in comoving 
units. Normally, we show explicit dependence on h (which de- 
notes i/o/100 kms" 1 Mpc" 1 ); but for absolute magnitudes we 
suppress this dependence, so that Mb denotes Mb — 5 log 10 h. 



2. Data base: The COMBO-17 Survey 

To date, COMBO-17 (Classifying Objects with Medium Band 
Observations in 17 filters) has surveyed three disjoint ~ 
31' x 30' southe rn equatorial fields (for their coordinates see 
IWolf et alJ2003l) to deep limits in 5 broad and 12 medium pass- 
bands, covering wavelengths from 400 to 930 nm. A detailed 
de scription of the su rvey along with filter curves can be found 
in IWolf et alJ J2004I) . All observations were carried out using 
the Wide Field Imager at the MPG/ESO 2.2 m-telescope on La 
Silla, Chile. 

In each filter, typically 10 to 20 individual exposures 
were taken (up to 50 for ultradeep i?-band images totalling 



20 ks with seeing ^ 0'.'8). Galaxies wer e detected on the 
deep .B-band images by using SExtractor JBertin & ArnoutsI 
Il99d) . The spectral energy distributions (SEDs) for i?-band de- 
tected objects were measured by performing seeing-adaptive, 
weighted-aperture photometry in all 17 frames at the position 
of the i?-band detected object. All magnitudes are quoted with 
a Vega zero point. 

Using the 17-band photometry, objects are classified us- 
ing a scheme based on template spectral energy distributions 
(Wo lf et alJ l2001blal) . The classification algorithm basically 
compares the observed colours of each object with a colour 
library of known objects. This colour library is assembled 
from observed and model spectra by synthetic photometry per- 
formed using an accurate representation of the instrumental 
characteristics of COMBO-17. Fo r galaxy classification, we 
use P EGASE model spectra (see iFioc & Rocca-Vo lmerange 
1997 for an earlier version of the model). The template spec- 
tra are a two-dimensional age/reddening sequence, in which a 
fixed exponential star formation timescale t = 1 Gyr is as- 
sumed, ages vary between 50Myr and 15 Gyr, and the redden- 
ing E(B — V) can be as large as 0.5 mag, adopting a Small 
Magellanic Cloud Bar extinction curve. Note that we do not 
apply any morphological star/galaxy separation or use other 
criteria. 

Using a minimum variance estimator, each object is as- 
signed a redshift (if it is not classified as a star). The redshift 
errors in this process depend on magnitude and type of the ob- 
ject, and for galaxies can be approximated by 

1/2 

(1) 



0.007 i + io°-8(«- 21 - 6 ) 



(1 + *) 

The galaxy redshift estimate quality has been tested by com- 
pari son with spectros copic redshifts for almost 1000 objects 
(see Wol f et alJl2004 . At bright limits R < 20, the redshifts 
are accurate to a z / (1 + z) ~ 0.01, and the error is dominated 
by mismatches between template and real galaxy spectra. This 
error can contain a systematic component that is dictated by 
the exact filter placement, but these 'redshift focusing' effects 
are of the order of magnitude of the random redshift errors for 
z < 1 and are unimportant for the current analysis. At the me- 
dian apparent magnitude R ~ 23, <x z /(l + z) ~ 0.02. For 
the faintest galaxies, the redshift accuracy approaches those 
achievable using traditional broadband photometric surveys, 
a z j (1 + z) ^ 0.05. We thus restricted our analysis to galaxies 
with I < 23. 

Fig. [2 shows the redshift distribution of the 22 310 
COMBO-17 galaxies between z = 0.2 and z = 1.2 (with 
I < 23 and M B < -18). The peak at z = 0.733 in Fig. [T] is 
due to a real structure in the Chandra Deep Field South, which 
has been spectroscopically confirmed JGilli et al.l2003l) . In or- 
der to define a volume limited sample, we restrict our analysis 
to the redshift range 0.4 < z < 0.8 and galaxies brighter than 
Mb — —18, which leaves us with 10 360 galaxies for the anal- 
ysis. Note that £?-band luminosities can be determined directly 
without any if-correction uncertainty, based on the photometry 
in our 17 filters between 400 and 930 nm and an interpolation 
of the corresponding template spectra. We do not apply any 
evolutionary corrections. 
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Fig. 1. The redshift distribution for galaxies with / < 23 and 
M B < -18 in the three COMBO-17 fields. The dotted lines in- 
dicate the redshift range over which we define a volume-limited 
sample for the clustering analysis. 

Table 1. The number of red sequence and blue cloud galaxies 
(by the definition of equations|3|5]and0, with I < 23, M B < 
-18 and 0.4 < z < 0.8 per COMBO-17 field. 



COMBO-17 field 


N Ic d 




CDFS 


752 


2782 


A901 


836 


2563 


Sll 


816 


2611 



The distribution of the redshift errors for all galaxies in our 
subsample with / < 23, M B < -18 and 0.4 < z < 0.8 is 
shown in Fig. [2] We use the prescription of B ell et all |2004) 
to separate galaxies into the red-sequence component and the 
remaining blue cloud component: 

Red sequence : (U — V) > (U — V)n m (2) 
Blue cloud : (U - V) < (U - V) hm (3) 
(U - V) iim = 1.25 - 0.4z - 0.08(Afy - 5 log 10 h + 20) , (4) 

where z denotes the redshift of each single galaxy. 

Note that the cut that separates the red sequence galaxies 
from the blue ones depends on both redshift and absolute V 
magnitude. This yields 2404 and 7956 galaxies in the red and 
blue subsamples; the former tend to have more accurate red- 
shifts, as shown in Fig. [2] This is not because the classification 
scheme works better for the red galaxies, but because they are 
on average brighter than the blue ones. 

Tabled shows the number of red sequence and blue cloud 
galaxies per COMBO-17 field. 

3. The halo model of galaxy clustering 

In discussing the results of clustering analyses of the COMBO- 
17 data, we will make frequent comparisons with theoretical 
predictions. We therefore now summarise the framework used 
to carry out this modell ing. This g oes back to the paradigm in- 
troduced bv lWhite & Reesl |l978): galaxies form through the 




0.02 0.04 0.06 0.08 0.1 0.12 

<7 Z /(1 + Z) 



Fig. 2. Distribution of relative redshift errors of galaxies with 
I < 23, M B < -18 and 0.4 < z < 0.8 (solid line). Also 
shown are the redshift error distributions for red sequence and 
blue cloud galaxies. 

cooling of baryonic material in virialized haloes of dark mat- 
ter. The mass function and density profiles of these haloes can 
be expressed in terms of simple fitting formulae derived from 
N-body simulations, and the large-scale clustering of haloes 
can be derived analytically for Gaussian density fields. This 
concentration on dark-matter haloes gives concrete form to 
earlier work on the clustering stati stics generated by dist ribu- 
tions of extended clumps (see e.g. iNevman & Scottlll952l and 
IScherrer & Bertsching erll99ll) . 

With an accurate description of dark-matter clustering to 
hand, the sta ge was set for an extension to galaxies via the 
'halo model' iMa & Frvll2000HSeliakll2000HPeacock & Smithl 
2000; iCoorav & S heth 2002). In this picture, the key remain- 
ing uncertainty is the way in which galaxies occupy the dark 
matter haloes; this can be regarded as an unknown function, to 
be probed experimentally. The halo model then allows a simple 
and direct understanding of many features of galaxy clustering, 
and how the clustering of galaxies differs from that of the mass. 

In this approach, the density field is a superposition of dark- 
matter haloes, with small-scale clustering arising from neigh- 
bours in the same halo. The corresponding real-space correla- 
tion function can be written as a combination of two parts: 

£r = Clin + 6ialo , (5) 

the first term representing the clustering of the dark matter 
haloes, and the second correlations from within a single halo. 
Large-scale halo correlations depend on mass, and the linear 
bias parameter for a given class of haloes, b(M), depends 
on the rare ness of the fluctuation and the rms of the under- 
lying field dKaiserlll984t ICole & Kaiserlll989i: IMo & White! 
1996; Sheth & Tormen 1999), usually measured in spheres of 
8/i _1 Mpc and termed as- The mass profile of the haloes is 
known from simulation s, and may be assumed to follow ei- 
ther a n NF W profile dNavarro et alJI 19961 Il997t iNavarro et ail 
|2004 . or a lMoore et al.Nl999l) profile. 

The key feature that allows bias to be included is to encode 
all the complications of galaxy formation via an halo occupa- 
tion number: the number of galaxies found above some lumi- 



4 



S. Phleps, J. A. Peacock et al.: Galaxy clustering from COMBO-17 



nosity threshold in a virialized halo of a given mass M. A sim- 
ple but instructive model for this halo occupation distribution 
(HOD) is 



N{M) 



= 1° 

\ (M/M c ) a 



(M < M c ) 
(M > M c ) 



(6) 



This is closely related to the mass-dependent weight introduced 
bv ljing et alJ Jl998l) . A model in which light traces mass ex- 
actly would have M c — ► and a = 1. The galaxies are as- 
sumed to be split into a central galaxy plus some number of 
satellite galaxies, which follow the mass distribution in the 
halo. It is necessary to make an assumption about the statis- 
tics of the HOD - in particular whether N is a causal function 
of M or whether it obeys a Poisson distribution. It is known 
that sensible results require sub-Poisson behaviour, and we as- 
sume the extreme limit in which N is perfectly determined by 
M, Putting all these ingredients together, the galaxy correlation 
function can be calculated analytically. 

Since its initial development, the halo model has been 
applied successfully to the interpretation of the correla- 
tion function of ga laxies in the local universe, notably by 
Zeha vi etalJ §004), who detected an inflection in the cor- 
relation function of SDSS red galaxies, interpreting it as 
indicating the transition regime between clustering domi- 
nated by 1-halo and 2-halo terms. The occupation mod el has 
been elaborated quite significantly (e.g . Abazaiia n et al. 2005; 
Zheng et all2005UKravtsov et alJ2004t iTasitsiomi et all2 004 ; 
Zent ner et alJl2005l) . including up to three parameters for the 
occupation distribution, plus the inclusion of some nonlinear 
evolution of the power spectrum in the £u n term. These sophis- 
tications can improve the detailed fit to correlation-function 
data from simulations, but they do run somewhat counter to the 
original heuristic spirit of the model. In this work, we shall re- 
tain the original m ethod of calculation, as described in detail by 
IPeacock & Smith! J2000l) . together with the simple power-law 
occupation model. This approach seems justifiable in a first ex- 
ploration of intermediate-redshift clustering, and the main fea- 
tures of interest are in any case relatively robust. 

The form of any transition-regime feature in the correlation 
function depends mainly on the mean halo mass occupied by 
the galaxies, - i.e. the average value of M, weighted by N{M) 
- and is insensitive to the details of their distribution within 
the halo. This typical halo mass is determined by our two- 
parameter model for N(M) plus the halo mass function. We 
shall use the additional constraint of the observed number den- 
sity of galaxies under study, so that there remains a single free 
parameter in the model - which we take to be the occupation 
slope a. For a given value of a, the number density determines 
the cutoff mass and hence the average halo mass. Different 
models of the HOD can of course be used; Zehavi et al. (2004) 
take N = 1 above M c to represent central galaxies, plus a 
power-law N(M) representing satellites, which commences at 
a mass of approximately 20M C with a slope of a ~ 1. In prac- 
tice, this model gives results similar to the single power law 
with a ~ 0.6, showing that the detailed shape of the HOD is 
hard to measure given only correlation-function data. The typ- 
ical number-weighted halo mass is a more robust quantity, and 



this is probably the best way to compare different HOD mod- 
els. 

We will normally assume a standard flat cosmology with 
n m = 0.25, fa = 0.045 and h = 0.73 where ff n = 100/ikm 
s" 1 Mpc" 1 ( Soe rgel et alJl2006t ICole et alJl2005h . The most 
uncertain cosmological parameter is the normalization of the 
power spectrum, cr§. We will often assume a standard value of 
(78 = 0.9, but it is also of interest to leave both erg and the 
power-law index a of the halo occupation number (equation|6j 
as free parameters, and determine them from a fit to the mea- 
sured galaxy correlation functions. 

4. Redshift space correlations 

4. 1. Projected correlations 

With a sample of the present size, only the two-point corre- 
lation function, £(r), can be determined accurately. Even this 
is not straightforward, because of the need to work in redshift 
space. While it is possible to measure angular positions on the 
sky with high precision, peculiar velocities as well as redshift 
errors distort the galaxy pattern along the line of sight, making 
£(r) appear anisotropic, and tending to reduce its amplitude. 

These problems can be overcome by splitting the separa- 
tion vector r of a pair of objects into components lying on the 
plane of sky, r p , and along the line of sight, n, and compute the 
correlation function £ (r p , ir) as a function of these two compo- 
nents. Projecting £(r p ,7r) onto the r p axis gives the projected 
func tion w(r„), which i s independent of any radial distortions 
(Dav is & Peeblesll 19831) . For small angles r 2 = r 2 + ir 2 . Thus 
the projected correlation function is defined as 



w(r p 



oc 
oo 



( 2 , 2\l/2 

(r p + 7T ) 1 



dn 



tkr){r 2 -r 2 p )- l l 2 rdr 



(7) 



Note that w(r p ) has dimensions of length. Since the correlation 
function converges rapidly to zero with increasing pair separa- 
tion, the integration limits do not have to be ±oo, but they have 
to be large enough to include all correlated pairs. As will be ex- 
plained in the following paragraph, this is a crucial point when 
the redshift errors are large. 

4.2. The effect of redshift errors on £(r p , ir) 

We now illustrate the effects of redshift errors, starting from a 
model for the true £(r p , ir). This is calculated using the halo 
model, as described above in section [51 together with a pre- 
scription for redshift-space distortions. Seliak (2001) showed 
how to include redshift-space distortions in the halo model, but 
it is also common to use the following simple model for the 
ratio of redshift-space power to real-space power: 



Ps/Pr = 



(1 + /3M 2 ) 2 
(l + fcV^/2) : 



(8) 



where o p is an effective pairwise velocity dispersion, /i is 
cosine of the angle between k and the line of sight, and 
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Fig. 3. The redshift-space correlation function £ (j 
lated using the halo model, plotted as a function of transverse 
(r p ) and radial (tt) pair separation. The data from the first quad- 
rant are repeated with reflection in both axes, in order to clarify 
deviations from isotropy. Upper panel: before convolution with 
the pairwise redshift error distribution inferred from the sample 
of COMBO-17 galaxies under consideration in this analysis, 
lower panel: after convolution. 



sequence of large-scale coherent motions. The former effect is 
not clearly visible owing to the large scale of the plot. 

The effect of redshift errors on this redshift-space correla- 
tion function is straightforward: it is a convolution in the radial 
direction. This reflects the fact that 1 + £ is a ratio of the ob- 
served and expected numbers of pairs of galaxies. The correla- 
tion function £(r p , tt) becomes distributed more broadly along 
the tt axis, but the total correlation signal is conserved - this is 
why w(r p ) is independent of redshift errors. In order to model 
this process, we need a model for the convolving function. This 
is simply deduced, because the classification scheme automat- 
ically retur ns an estimate of the rms redshift error for each 
galaxy (see IWolf et al.ll2004l) . Thus, given a pair of galaxies 
i, j, with redshift errors a Zi and a Zj the rms pairwise error is 
(Tp a i r . . = (a z . + al ) 1 / 2 . The signal from this pair is smeared 
by a Gaussian with this width, so the overall convolving func- 
tion is a sum of the Gaussians corresponding to all pairs: 



f(5z) = ^(Znap^y 



1/2 



exp [-(6z/a pah%n ) 2 /2] , (9) 



where N is the number of pairs After having transferred 
the pairwise redshift error distribution into comoving distances 
we can convolve £(r p , tt) with Eq. 

The effect of this convolution is shown in the second panel 
in Fig. |3] The redshift-space correlations are now heavily elon- 
gated in the radial direction, and some care is needed in extract- 
ing the projected correlation signal. 




(3 = fl^f/b (e.g. Ballin ger. Peacock & Heavenslll996l) . This 
approach has the advantage that the redshift-space correlation 
function can then be foun d analytically, apart from a radial con- 
volution (Hamilton 1992). We used this prescription, taking a p 
to be \/2 times the one-dimensional velocity dispersion cal- 
culated from the halo model. Fig. [3] shows a model £(r p ,7r), 
which contains two well-known expected anisotropies: the iso- 
correlation contours of £(r p ,w) are stretched along the tt di- 
rection at small separations, because of the effect of virialized 
velocity dispersions, and compressed at large scales as a con- 



Fig. 4. An example of the model-fitting procedure used to de- 
termine w(r p ): The grey line shows the correlation function 
£(r p , tt) of the red sequence galaxies as a function of tt at 
r p — 0.2 hr 1 Mpc. The dashed line shows the model corre- 
lation function in the absence of redshift errors. The dotted and 
solid lines show the smoothed data (the smoothing was done 
by applying a box filter of length 40 hr 1 Mpc along the 7r-axis 
using the running average), and the model after the convolution 
with the expected error distribution, with the solid part indicat- 
ing where the integration is carried out. 
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4.3. Estimation of projected correlations 

The simplest strategy for carrying out the projection needed 
in order to deduce w(r p ) would be to integrate £,(r p ,ir) over 
a very large radial range. Fig. [3] suggests that a maximum n 
value of 150 to 200/i -1 Mpc would be required to capture 
all the signal. The problem with this strategy is that the ran- 
dom noise in £(r p , ir) is independent of it at a given r p (be- 
cause the expected pair counts have a cylindrical dependence 
oc r p dr p dir). Thus, integration to ir = 200/i _1 Mpc would 
yield a random error in w(r p ) that is \/2 times larger than in- 
tegration to 7r = 100 h~ l Mpc - but the lower limit systemati- 
cally misses part of the signal. 

We have developed a strategy for solving this problem, 
which depends weakly on some prior knowledge of the likely 
form of the true clustering signal (after error convolution). A 
model for £(r p , it) defines how the real-space signal w(r p ) is 
spread out in 7r; we are only concerned with the shape of this 
function, which is dominated by the convolution with the red- 
shift error distribution. Given this probability distribution in the 
7r direction, the amplitude of w(r p ) can be estimated by fitting 
a scaled version of our model £(r p , it) to the data at the r p value 
of interest. 

In practice, our model (,(r p ,ir) will not be exact, and we 
considered the following compromise procedure for estimat- 
ing w(r p ) so that the result is robust. For each r p value, we 
fit the amplitude of the (convolved) model £(r p — const, 7r) 
to the data. We then integrate the data for £(r p ,7r) out to 
7r = 100 h^ 1 Mpc, from which point on we integrate the con- 
volved model out to infinity (see Fig. 0). This combines the 
exact measurement of w(r p ) within 7r max = 100 /i -1 Mpc 
with an estimate of the missing signal at larger ir. Since this 
correction is typically 20% of the overall signal, we do not 
need to estimate it very accurately. In practice, the results from 
the 2-stage procedure were very similar to the direct fitting 
method. This process is performed separately for the red and 
blue galaxy samples, using the appropriate pairwise error dis- 
tributions and the final best-fitting halo-model £(r p , it) for the 
unconvolved prediction. 

The width of the convolved model in Fig. [4] suggests that 
the redshift errors yielded by the object classification scheme, 
which we used for the calculation of the pairwise error distri- 
bution, may be slightly overestimated. In order to estimate the 
effect on the projected correlation function, we tried repeating 
the analysis with redshift errors scaled to 80% of the values 
given in the object catalogues. This scaling gives the best fit to 
the data; however, the resulting changes to w(r p ) were small 
compared to the random errors. 

4.4. Integral constraint 

The mean galaxy density is determined from the observed 
galaxy counts in each field, whic h does not ne cessarily repre- 
sent the the true density (Groth & Peebles 1977). The estimator 
will be on average biased low with respect to the true correla- 
tion by a constant X: 



where Wt (r p ) is the true projected correlation function, w m (r p ) 
the measurement. The integral constraint X is given by 



1 — J2 / w t (r p )d 2 Sxd 2 S 2 



(11) 



where S is the physical area corresponding to the solid angle of 
the field at the redshift under consideration. For the calculation 
of the integral constraint, we assume that the three dimensional 
correlation function £(r) is to first approximation a power law: 



(12) 



£(>") = - 

Jo 



Then the evaluation of equationQyields 

w(r p ) = Crlrl' 1 , 



(13) 



where C is a numerical factor, which depends only on the slope 



7: 



c = vtF 



r((7-i)/2) 



(14) 



r(7/2) ' 

If the true correlation function is given by equation Jl 3I >. the 
measurement yields 



w m [r v 



or r p 



CrZ[r l -"i-l/{CrZ)] 



(15) 
(16) 



The true amplitude CrJ is not known, but XjiCr^) can be 
estimated by performing a Monte Carlo integration (where we 
use the mean of the pair counts (RR) at a projected distance r p 
of the four fields): 



1 _ g [(RR) ■ r£T] 



C 



(17) 



w m {r p ) = w t (r p ) 



(10) 



The true value of CYq can be estimated by fitting equation dl6l 
to the data, taking the value of X/(CYq) from equation d!7i . 
This value, multiplied by the fitted amplitude CYq, yields the 
integral constraint X. The measurement can then be corrected 
for the integral constraint by adding X to w m (r p ). 

This method yields estimates of X = 0.14 h^ 1 Mpc for 
the red galaxies and X — 0.33 ft. -1 Mpc for the blue galax- 
ies. These values are negligible in comparison with the ob- 
served data for r p 20 h^ 1 Mpc, demonstrating that the fields 
are large enough to deliver a fair sample. As a cross-check, 
note that we expect X = 200er 2 (since we integrated over 
An — 200 ft -1 Mpc), where 1 is the fractional variance in 
galaxy numbers between different realizations of our survey. 
With three fields, a should be y3 times smaller than the field- 
to-field rms variation, so our figures for X suggest 4.5% and 7% 
expected scatter in the numbers of galaxies per field for red and 
blue galaxies respectively. This agrees well with the numbers 
in Table 1 . 

4.5. Error analysis 

Finally, there is the crucial issue of setting realistic error bars 
on our correlation estimates. The three COMBO-17 fields mea- 
sure ~ 31' x 30' each and are thus large enough to carry out 
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Fig.5. £(r p , tt) of all COMBO-17 galaxies with 0.4 < z < 0.8, 
/sis < 23 and Mb < —18. Again the data from the first 
quadrant are repeated with reflection in both axes. In the trans- 
verse direction the pair separations accessible for the analysis 
are limited by the COMBO-17 field of view. 

a jack-knife analysis. We divide each field into four quadrants, 
and then calculate the correlation function w(r p ) (including the 
integral constraint) for twelve realisations of the data, each time 
omitting one of the quadrants. The variance in w is then given 
approximately by 



N-l 

N 



E 

= 1,JV 



>(>>)) 



Wi(r v 



where N = 12 is th e number of realisations of the data (e.g. 
IScranton et al.hOOA. 

In order to check for cross-correlations between the data 
points, we can extend the jack-knife method in the obvious way 
to estimate the covariance between different bins, cr 2 , . The nat- 
ural way to express this is as a correlation coefficient matrix: 
nj = afJaiCTj. Results in this form are presented below. 



5. The clustering of the COMBO-17 galaxies 

5. 1. Results 

We calculated £(r p , tt) for all COMBO-17 galaxies in the red- 
shift range 0.4 < z < 0.8 with /-band magnitudes / < 23 and 
absolute restframe B ban d luminosities Mr < — 18. We used 
the estimator invented bv lLandv & Szalavl lll993h . An angular 
mask for the survey was derived by censoring the surround- 
ings of bright stars in the fields. The same mask was applied to 
a random catalogue consisting of 30 000 randomly distributed 



galaxies, each of which was assigned a redshift taken randomly 
from the real data, where the three fields were put together 
in order to smooth the redshift distribution. Using a smoothed 
form of the empirical redshift distribution did not yield a sig- 
nificant change in the results. 

The resulting £(r p , tt) is shown in Fig. [5] The field of view 
of the COMBO-17 fields limits the pair separations accessible 
for the analysis, so in the transverse direction there is of course 
no signal at separations larger than the physical distance corre- 
sponding to the diagonal diameter of the fields. 

For each object, we have an estimate of the redshift and the 
restframe colours and luminosities; it is therefore possible to 
divide the sample into two distinct colour classes as described 
earlier. For both samples we calculated w(r p ) as described in 
section[4l correcting for the integral constraint I, and the influ- 
ence of the redshift errors. These results are shown in Fig. [6] 

5.2. Fitting the halo model 

Fig. [6] also shows predictions from the halo model, varying 
the single occupation-number parameter a, and choosing the 
cutoff M c so as to match the observed comoving densities of 
0.004ft, 3 Mpc" 3 (red) and 0.012/i 3 Mpc" 3 (blue). It is appar- 
ent that there is greater sensitivity to a at small separations, 
and that once a is fixed from the data there, there is little 
freedom at large separations, where the data and the model 
match satisfyingly well. The preferred values are approxi- 
mately a = 0.5 for the red population and a — 0.2 for blue 
galaxies. These figures correspond to cutoff masses of respec- 
tively M c = 10 12 - 15 h- 1 M Q and M c = lO 11 - 50 ^ 1 M Q . As 
discussed earlier, a more meaningful way of casting these num- 
bers may be to apply the HOD model to the halo mass function, 
to calculate the effective halo mass, weighting by galaxy num- 
ber. These figures come out as M e g = 10 13 21 /i~ 1 Af© and 
M eS = 10 12 - 52 /i _1 M Q respectively. 

Fig. [6] also shows a magnified view, with the measured 
correlation functions and the corresponding best-fitting models 

< 



(18) both divided by a power-law fit (fitted in the range log 



10 ' p 

1.1), the slope and amplitudes of which are given in Table [2] 
The data points do not scatter arbitrarily around the power-law 
fit, but show systematic deviations. For the red galaxies, there 
is a marked dip around r p ~ 1.5 h^ 1 Mpc; the blue galaxies 
are closer to a power law, but with a relatively abrupt step at 
r p ~ 0.2 h^ 1 Mpc. Both these features are impressively well 
accounted for by the halo model predictions, especially when 
it is considered that there is only one free parameter. 

It is interes ting to compare our r esults with those of the 
VVDS project jLe Fevre et aTll2005l) . They give results to a 
similar depth for two fields, although not divided by colour, 
with a total of 7155 redshifts over 0.61 deg 2 . Their redshift 
bins are not identical, but they quote tq = h~ 
and 7 = 1.71±£^ at (z) = 0.6 and r = 4.55^ , !; 



,+0.18 
-0.11 

and 7 = 1.48±g ?| at (z) 
CDFS, which is one of our fields, and we have checked that our 
figure for this field alone agrees well with the VVDS, as it does 
for our other fields. The VVDS 2 h field thus gives a somewhat 
lower clustering strength; this may be because the VVDS sam- 



+0.53 ^-l Mpc 

1 25 ft- 1 Mpc 
0.7. The latter figure is from the 
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Table 2. The amplitudes, slopes and resulting correlation 
lengths r of the best fitting power law (w(r p ) = Ar p , fitted 
for log 10 r p < 1.1), for the total sample, and divided into red 
sequence and blue cloud galaxies. 









53.98±i;gg 


-0.57 ±0.02 


4.n±8:S 




Srcd 


r r cd 


86.28±US 


-0.94 ± 0.03 


5.39±8;|g 


■Ablue 


<5bluc 


r blue 


37.09±|;|f 


-0.45 ± 0.03 


3.64±S!;i 



pie in that field is about 0.5 mag. deeper than the one studied 
here, plus the fact that the VVDS analysis has no lower limit in 
luminosity. The ratio in ro between the VVDS 2 h field and our 
overall result is 1.75 ± 0.21, so a factor 1.3 from luminosity- 
dependent clustering would be required in order to make the 
results statistically consistent. 

We now want to quantify the agreement between model and 
data in more detail, using the jack-knife error estimates. This 
would be straightforward if the covariance matrix were diago- 
nal, and if the model was exact. In practice, there is some de- 
gree of correlation in the data, and the simple halo model used 
here may be expected to have some systematic deviations with 
respect to ideal data. The correlation matrices in Fig.[7]indeed 
show a strong correlation between the large-scale data points 
(large matrix indices, lower right corner). This is due to the 
integral constraint, which is included in the calculation of the 
jack-knife realisations: all large-scale data points are offset by 
the same amount and thus become correlated. This issue is not 
too serious, since the errors in this regime are in any case large. 
We therefore ignore the large-scale points at r p > 10 hr 1 Mpc 
and treat the remaining data as independent. Even so, our ap- 
proximate model is not guaranteed to deliver a perfect fit. In 
order to achieve a formally acceptable value of \ 2 for the fit, 
we added in quadrature an error of 5% to the errors on w for 
red galaxies. When fitting the z = data, as discussed below, a 
covariance matrix was not available, and the formal errors are 
in any case small. In this case, we therefore took what is ef- 
fectively a least-square approach, which required an effective 
error of 10% in both blue and red galaxies. 

In performing the fitting, it is interesting to consider vari- 
ations in both the power-law index a of the HOD (equation 
|6jl, and in the normalization erg. We emphasise that erg 1S the 
zero-redshift value, which is connected to the degree of inho- 
mogeneity at z = 0.6 by the growth factor predicted by the cos- 
mological model (a change in linear density contrast by a factor 
1.32). Fig. [8] shows the likelihood contours for these two free 
parameters a and erg. The preferred values and marginalized 
rms errors are a — 0.56±0.03 and erg = 0.84±0.08 for the red 
sequence galaxies and a = 0.16 ± 0.03 and erg = 1.19 ± 0.09 
for the blue population. 



These independent estimates of erg from the red and blue 
populations are both close to our default value of 0.9. The 
agreement is not perfect, and the difference in a% is formally 
3cr, but the mean value of <7g = 1.02 ± 0.17 is certainly plausi- 
ble. Given the simplicity of the HOD model, this is a satisfying 
result. 

5.3. Robustness of the results 

In view of the tension between the normalization inferred from 
the red and blue galaxies, and in the light of the preference 
for a low normalization of as — 0.75 from the 3 -year WMAP 
data ( Sper gel et alJl 2006). it is important to discuss the extent to 
which our result is rendered uncertain by simplifying assump- 
tions in the modelling. 

The halo model is an idealized approximation in many 
ways, and one issue in particular has generated considerable 
discussion recently. It has been shown that the clustering of 
dark mat ter haloes is dependent on the halo formation time 
dSheth fc Tormenl l2004alri lOao et al I EoTfl IWechsler et al. I 
2005| lHarker et alJl2006l: iReed et alJEoOri) . and the question 
is what impact this has on halo model calculations. To some 
extent, the effect is already included in the halo-model for- 
malism: when the cosmic density field is smoothed on a given 
mass scale, the clustering of peaks in the smoothed field is well 
know n to increase with peak height, i.e. with formation red- 
shift jKaiseJ fT984). More massive systems are more strongly 
clustered for the same reason: only the rarest peaks exceed 
the threshold for collapse when the variance in filtered den- 
sity is low. When we use the standard expression for bias as a 
function of mass, this averages over all systems that have col- 
lapsed by the present: the dependence of clustering on forma- 
tion time will thus have no effect on predicted galaxy prop- 
erties if the occupation numbers are purely a function of halo 
mass. However, it seems reasonable that the occupation num- 
ber for a given mass will in fact depend to some extent on col- 
lapse redshift (e.g. more red galaxies in a halo that collapses 
early). At a minimum, the age-clustering effect will then con- 
tribute to a stochastic aspect of the occupation number, so that 
there is some scatter in N at a given Al. More seriously, it can 
also bias the mean clustering compared to all haloes of that 
mass. The influence of these effects on halo-model predictions 
remains to be explored, and this task is beyond the scope of 
th e present paper . Some work a long these lines has be e n don e 
bv lZentneretai] (|2005) and bv lCroton. Gao & White! ll2006l) . 
where the halo contents in a simulation are scrambled between 
all haloes of the same mass, thus de s troyin g any correlations 
with collap se redshift. |Zentner et al. (2005)" did this for sub- 
haloes, but ICroton. Gao & W hite (2006) considered the case 
of most direct interest, which is semianalytic galaxy popula- 
tions. They do detect systematic shifts in correlation ampli- 
tude, but for the luminosities of interest here these are no larger 
than 10%. Such shifts are not important in comparison with the 
COMBO-17 measuring errors, but this is clearly an issue that 
should be looked at in more detail, especially as the accuracy 
in measuring high-redshift clustering improves. 
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Table 3. The inferred values of a s from COMBO-17 red and 
blue galaxies, showing the effect of variations in some of the as- 
sumptions in the halo-model calculation. We also give the cor- 
responding values of a, the power-law slope in the occumap- 
tion model, plus formal values of \ 2 for the best-fitting models 
(on 20 data points). 



We use the flux-limited 
and the -19 > 



Halo Model 


red 
°8 


blue 
OS 


Ctrcd 


Qbluc 


Xrcd 


2 

Xbluc 


Standard 


0.84 


1.19 


0.56 


0.16 


26.2 


9.7 


Sim = 0.2 


0.80 


1.07 


0.52 


0.16 


38.0 


16.1 


h = 0.6 


0.90 


1.15 


0.51 


0.14 


33.0 


10.4 


1 + 5 V = 300 


0.90 


1.02 


0.55 


0.18 


43.1 


9.6 


Zheng N{M) 


0.85 


0.98 


0.96 


0.56 


46.4 


9.4 



Other degrees of freedom in the halo model are more eas- 
ily investigated, and we summarise some tests here, the results 
of which are presented in Table [3] We show the impact on the 
values of erg inferred from red and blue galaxies by (a) vary- 
ing cosmological parameters; (b) varying parameters internal 
to the halo model; (c) varying the occupation number prescrip- 
tion. In the first category, we see that the Hubble parameter has 
very little effect, but that <7g increases with il m very roughly as 
rj° rl 3 to £!° n D , with a larger sensitivity for the blue galaxies. In 
the second category, we considered altering the assumed den- 
sity contrast for a virialized halo from the usual figure of 200 
to a slightly larger number, which is sometimes assumed in a 
low-density model. Finally, we modify our simple power-law 
N(M) to something that resembles more closely the prescrip- 
tion used by e.g. lZheng et alJ J2005I) : N = 1 between M c and 
10Af c , rising as M a for smaller M (we considered a < 1). 
We also include the formal \ 2 values for some of these alter- 
natives. It appears that the standard model provides the best fit, 
especially to the red galaxies. However, bearing in mind the 
simplified nature of the halo model, these differences should 
not be given high weight; it is more interesting to concentrate 
on the robustness of the best-fitting parameters. 

The overall conclusion of these tests is that plausible varia- 
tions of some of the degrees of freedom in the modelling can al- 
ter erg by 10 to 20%, and in a way that changes the consistency 
between red and blue results by a similar amount. We there- 
fore conclude that the modelling is working as well as could 
have been expected, and that there is no need to be concerned 
by either the internal red-blue tension or by the 1 .6a discrep- 
ancy with the WMAP erg . We now turn to the comparison with 
z = 0; some of the systematics in the analysis should be com- 
mon to all redshifts, so we should hope for a good level of con- 
sistency between measurements based on data from different 
epochs. 

5.4. Comparison with local clustering 

As a local comparison, we considered w(r p ) for a combined 
set of red and blue galaxie s taken from both Sloan Digita l Sky 
Surve y (SDSS; lYork et alj EoOO) and 2dFGRS JColless et al l 
12001 1) data. The sample has been divided into red and blue by 
either the bimodality of the rest-frame colour distribution, or 
spectral type, in a way that should compare reasonably well 



with the COMBO-17 classification. 
2dFGRS results of lHawkins et al l d2003ah . 
M r > -20 results of lZehavi et alJ ll2005l) . which have closely 
comparable amplitudes. Fig.|9]shows a comparison of the local 
and the high redshift sample. What is apparent here is that the 
main difference between the COMBO-17 sample at z = 0.6 
and the local data is in the shape of the correlation function, 
with almost identical amplitudes at small scales, but a differ- 
ence of nearly a factor two for r p > \ hr Y Mpc. 

As shown in Fig.|9] the halo model is capable of accounting 
for these differences. Fig.[H)|repeats the model-fitting exercise 
for the z = data, and the preferred values and marginal- 
ized rms errors are a = 0.49 ± 0.02 and er 8 = 1.03 ± 0.07 
for the red sequence galaxies and a = 0.23 ± 0.02 and 
as = 1.00 ± 0.05 for the blue population. These values of 
a are slightly smaller than those obtained from COMBO- 
17; they correspond to cutoff masses of respectively M c = 
10 i2.2o ft -i MQ and Mc = iq11-50 /i -i Mqi or effective halo 

masses M eff = lO 13 - 50 ^ 1 ^/© and M cS = lO 1280 /^ 1 M 
respectively. Compared to the z — 0.6 results, M c has not 
changed much, but AI c g has increased by about a factor 2. 
This makes sense in terms of hierarchical growth: the minimum 
dark-matter mass needed to assemble a galaxy-sized amount of 
baryons should be invariant, but such haloes inevitably merge 
into larger systems as time progresses. 

The mean value of cr 8 — 1.01 ± 0.04 from the z = data 
agrees very well with a s = 1.02 ± 0.17 from COMBO-17. 
This agreement assumes hierarchical growth in the halo mass 
function between z = 0.6 and the present, without which the 
inferred values of erg would have been expected to differ by a 
factor 1.3. 



6. Summary & conclusions 

Using a sample of 10 360 galaxies with photometric redshifts 
from the COMBO-17 survey, we have investigated in some de- 
tail the shape of the correlation function at redshift (z) ~ 0.6. 
We have shown for the first time that the two-point correlation 
functions of both red sequence galaxies and blue cloud galaxies 
at this redshift display deviations from a power law, analogous 
to the deviations seen at low redshift ( lHawkins et al.l l2003a: 
IZehavi et alJ2004 . 

We have compared these observations to the predictions of 
a simple halo model, and find a good fit. It appears that the 
COMBO-17 data allow us to identify the point of transition be- 
tween 1-h alo clustering and 2 -halo clustering, as was done at 
z = by Zeh avi et alJ J2004). The implication is that the red 
and blue galaxies at (z) — 0.6 inhabit haloes of typical effec- 
tive mass M cS = 10 13 - 2 h _1 M Q and M eS = 10 12 - 5 h- 1 M Q 
respectively. 

We have also allowed the zero-redshift normalization of the 
power spectrum, og, to be a free parameter in this analysis. 
Impressively, both red and blue subsets imply a consistent lo- 
cal normalization of the power spectrum: as — 1- This figure is 
close to the value inferred by independent means using CMB 
and gravitational lensing (e.g. Refregier 2003), and certainly 
within the tolerance expected from the inevitable systematics 
associated with the simple modelling that we have used in this 
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first analysis of galaxy clustering at intermediate redshifts. This 
consistency is obtained by assuming that the dark halo mass 
function grows in a standard hierarchical fashion so that the 
normalization at z = 0.6 is approximately 30% lower than to- 
day. Our results amount to a verification that growth of this 
order has occurred. 

We intend to expand this work to higher redshifts using 
COMBO-17+4, the NIR extension of COMBO-17, for which 
observations are currently being carried out using the 2k x 2k 
Omega2000 camera at the 3.5m-telescope on Calar Alto, Spain. 
Combining the existing optical data base from COMBO-17 
with NIR observations in one broad and three medium band 
filters (covering the wavelength range from 1040 to 1650 nm), 
we expect to obtain ~ 4200 galaxy redshifts with an accuracy 
of a z j (1 + z) = 0.02 up to z — 2. This much longer baseline 
in cosmic time will allow us to observe much larger evolution 
of halo masses, testing the idea of hierarchical growth back to 
a time close to the formation of luminous galaxies. 

Acknowledgements. S. Phleps acknowledges financial support by the 
SISCO Network provided through the European Community's Human 
Potential Programme under contract HPRN-CT-2002-00316. JAP was 
supported by a PPARC Senior Research Fellowship. CW was sup- 
ported by a PPARC Advanced Fellowship. We thank Luigi Guzzo for 
helpful comments. We thank the anonymous referee for helpful com- 
ments and suggestions. 



References 

Abazajian, K., Zheng, Z., Zehavi, I., et al. 2005, ApJ, 625, 613 
Ballinger, W.E., Peacock, J.A. & Heavens, A.F., 1996, 

MNRAS, 282, 877 
Bell, E. E, Wolf, C, Meisenheimer, K., et al. 2004, ApJ, 608, 

752 

Benson, A. J., Cole, S., Frenk, C. S., Baugh, C. M., & Lacey, 

C. G. 2000, MNRAS, 311, 793 
Bertin, E. & Arnouts, S. 1996, A&AS, 117, 393 
Cen, R. & Ostriker, J. P. 2000, ApJ, 538, 83 
Cole, S. et al. (the 2dFGRS Team) 2005, MNRAS, 362, 505 
Cole, S. & Kaiser, N. 1989, MNRAS, 237, 1127 
Cole, S., Lacey, C. G., Baugh, C. M., & Frenk, C. S. 2000, 

MNRAS, 319, 168 
Colin, P., Klypin, A. A., Kravtsov, A. V., & Khokhlov, A. M. 

1999, ApJ, 523, 32 
Colless, M., Dalton, G, Maddox, S., et al. (the 2dFGRS Team) 

2001, MNRAS, 328, 1039 

Cooray, A. & Sheth, R. 2002, Phys. Rep., 3 72, 1 

Croton, D.J., Gao, L., White, S.D.M. 2006, astro-ph/0605636 

Davis, M. & Geller, M. J. 1976, ApJ, 208, 13 

Davis, M. & Peebles, P. J. E. 1983, ApJ, 267, 465 

Fioc, M. & Rocca-Volmerange, B. 1997, A&A, 326, 950 

Gao, L., Springel, V. & White, S. D. M. 2005, MNRAS363, 

L66 

Gilli, R., Cimatti, A., Daddi, E., et al. 2003, ApJ, 592, 721 
Groth, E. J. & Peebles, P. J. E. 1977, ApJ, 217, 385 
Hamilton, A.J.S., 1992, ApJ, 385, L5 

Harker, G, Cole, S., Helly, J., Frenk, C. & Jenkins, A. 2006, 
MNRAS, 367, 1039 



Hawkins, E., Maddox, S., Cole, S., et al. (the 2dFGRS Team) 

2003, MNRAS, 346, 78 
Jenkins, A., Frenk, C. S., Pearce, F. R., et al. 1998, ApJ, 499, 

20 

Jing, Y. P. 1998, ApJ, 503, L9 

Jing, Y. P., Mo, H. J., & Boerner, G. 1998, ApJ, 494, 1 
Kaiser, N. 1984, ApJ, 284, L9 

Kauffmann, G, Colberg, J. M., Diaferio, A., & White, S. D. M. 

1999, MNRAS, 303, 188 
Kravtsov, A. V. and Klypin, A. A. 1999 ApJ, 520, 437 
Kravtsov, A. V., Berlind, A. A., Wechsler, R. H, Klypin, A. A., 

Gottlober, S., Allgood, B. & Primack, J. R. 2004 ApJ609, 35 
Landy, S. D. & Szalay, A. S. 1993, ApJ, 412, 64 
Le Fevre, O., et al. 2005, A&A, 439, 877 
Ma, C.-P. and Fry, J. N. 2000 ApJ, 543, 503 
Meisenheimer, K. & Roser, H.-J. 1986, in Use of CCD 

Detectors in Astronomy, Baluteau J. -P. and D'Odorico S., 

(eds.), 227 

Mo, H. J. & White, S. D. M. 1996, MNRAS, 282, 347 
Moore, B., Quinn, T, Governato, F, Stadel, J., & Lake, G. 

1999, MNRAS, 310, 1147 
Navarro, J. F, Frenk, C. S., & White, S. D. M. 1996, ApJ, 462, 

563 

— . 1997, ApJ, 490, 493 

Navarro, J. F, Hayashi, E., Power, C, et al. 2004, ApJ, 349, 
1039 

Neyman, J. and Scott 1952, ApJ, 116, 144 

Neyrinck, M. C, Hamilton, A. J. S. & Gnedin, N. Y. 2004 

MNRAS, 348, 1 
Norberg, P., Baugh, C. M., Hawkins, E., et al. 2002, MNRAS, 

332, 827 

Peacock, J. A. & Smith, R. E. 2000, MNRAS, 318, 1 144 
Pearce, F. R., Jenkins, A., Frenk, C. S., et al. 1999, ApJ, 521, 
L992 

Peebles, P. J. E. 1974, A&A, 32, 197 

Phleps, S. & Meisenheimer, K. 2003, A&A, 407, 855 

Reed, D. S., Governato, E, Quinn, T, Stadel, J., Lake, G. 2006, 

astro-ph/0602003 
Refregier, A. 2003, ARAA, 41, 645 
Roser, H. . & Meisenheimer, K. 1991, A&A, 252, 458 
Scherrer, R. J. & Bertschinger, E. 1991, ApJ, 381, 349 
Scranton, R. et al., 2002, ApJ, 579, 48 
Seljak, U. 2000, MNRAS, 318, 203 
Seljak, U. 2001, MNRAS, 325, 1359 
Sheth, R. K. & Tormen, G. 1999, MNRAS, 308, 1 19 
Sheth, R. K. & Tormen, G. 2004, MNRAS, 349, 1464 
Sheth, R. K. & Tormen, G. 2004, MNRAS, 350, 1385 
Somerville, R. S., Lemson, G, Sigad, Y, et al. 2001, MNRAS, 

320, 289 

Spergel, D. N., Verde, L., Peiris, H. V., et al. 2003, ApJS, 148, 
175 

Spergel, D. N., Bean, R., Dore, O., et al. 2006, 

astro-ph/0603449 
Tasitsiomi, A., Kravtsov, A. V., Wechsler, R. H. & Primack, 

J. R. ApJ, 614, 533 
Totsuji, H. & Kihara, T 1969, PASJ, 21, 221 
Wechsler, R. H., Zentner, A. R., Bullock, J. S., Kravtsov, A. V. 

2005, astro-ph/05 12416 



S. Phleps, J. A. Peacock et al.: Galaxy clustering from COMBO-17 



11 



Weinberg, D. H., Dave, R., Katz, N., & Hernquist, L. 2004, 
ApJ, 601, 1 

White, S. D. M. & Rees, M. J. 1978, MNRAS, 183, 341 
Wolf, C, Meisenheimer, K., Kleinheinrich, M., et al. 2004, 

A&A,421,913 „ 
Wolf, C, Meisenheimer, K., Rix, H.-W.; Borch, A., Dye, S., = 

Kleinheinrich, M. 2003, A&A, 401, 73 I 
Wolf, C, Meisenheimer, K., & Roser, H. . 2001a, A&A, 365, ^ 

660 £ 
Wolf, C, Meisenheimer, K., Roser, H. ., et al. 2001b, A&A, 

365,681 

York, D. G., Adelman, J., Anderson, J. E., et al. 2000, AJ, 120, 
1579 

Yoshikawa, K., Taruya, A., Jing, Y. P., & Suto, Y. 2001, ApJ, 
558, 520 

Zehavi, I., Weinberg, D. H., Zheng, Z., et al. 2004, ApJ, 608, 
16 

Zehavi, I., Zheng, Z., D.H. Weinberg, D., et al. 2005, ApJ, 630, 
1 

Zentner, A. R., Berlind, A. A., Bullock, J. S., Kravtsov, A. V. 

& Wechsler, R. H.2005 ApJ, 624, 505 1 
Zheng, Z., Berlind, A., Weinberg, D., et al. 2005, ApJ, 633, 791 2 




r p / h 1 Mpc 




Fig. 6. The projected correlation function w(r p ) for the 
COMBO-17 data, divided into red and blue galaxies. The lines 
show predictions from the halo model with a = 0.45, 0.5, 0.55 
for red sequence and a = 0.15, 0.2, 0.25 for blue cloud galax- 
ies, respectively, where increasing a corresponds to increas- 
ing clustering. The lower panels show an expanded view of the 
data, divided by the best-fitting power law (fitted in the range 
log 10 r p < 1.1). 
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Fig. 8. The likelihood contours for the two free parameters a 
and (T 8 . Upper panel: red galaxies; lower panel: blue galaxies. 
Contours are shown at the Gaussian equivalent of the following 
confidence levels: 68% 1-parameter, 68%, 95%, 99%, 99.9%, 
99.99% 2-parameter, i.e. Aln£ = 0.5, 1.1, 3.0, 4.6, 6.9, 9.2. 



Fig. 7. The correlation matrices for red (upper panel) and blue 
(lower panel) galaxies, for —1.4 < log 10 r p < 1.3. The appar- 
ent strong correlation between the large scale data points (lower 
right corner) is due to the integral constraint. 



S. Phleps, J. A. Peacock et al.: Galaxy clustering from COMBO-17 



13 




r p / h 1 Mpc 




r p / h 1 Mpc 



Fig. 9. The projected correlation function in the redshift bin 
0.4 < z < 0.8 (filled circles) in comparison with the local 
data (open circles). The upper panel shows red sequence galax- 
ies, the lower panel is for blue cloud galaxies. The lines show 
the prediction of the halo model, which are higher in the z = 
case. 




Fig. 10. The likelihood contours for the two free parameters a 
and as, now fitting the z = data. Upper panel: red galax- 
ies; lower panel: blue galaxies. Contours are shown at the 
Gaussian equivalent of the following confidence levels: 68% 
1-parameter, 68%, 95%, 99%, 99.9%, 99.99% 2-parameter, i.e. 
Aln£ = 0.5, 1.1,3.0,4.6,6.9,9.2. 



